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We consider the flow-driven translocation of single polymer chains through nanochannels. Using 
analytical calculations based on the de Gennes blob model and mesoscopic numerical simulations, 
we estimate the threshold flux for the translocation of chains of different number of monomers. The 
translocation of the chains is controlled by the competition between entropic and hydrodynamic 
effects, which set a critical penetration length for the chain before it can translocate through the 
channel. We demonstrate that the polymers show two different translocation regimes depending on 
how their length under confinement compares to the critical penetration length. For polymer chains 
longer than the threshold, the translocation process is insensitive to the number of monomers in the 
chain as predicted in Sakaue et ai, Euro. Phys. Lett., 72 83 (2005). However, for chains shorter 
than the critical length we show that the translocation process is strongly dependent on the length 
of the chain. We discuss the possible relevance of our results to biological transport. 



I. INTRODUCTION 

The passage of polymer chains through nanochannels 
is an ubiquitous process in nature. Biopolymers, such 
as DNA and RNA, have to cross a multitude of barriers 
to perform different biological functions, for example, in 
translocating through cellular membrane pores or when 
ejecting from viral capsids Considerable interest has 
arisen in the details of the translocation process due to 
a vast array of practical applications that include the 
potential sequencing of DNA chains [2|, the sorting of 
biopolymers using smart entropic traps Q or in sieving 
processes Q in the pharmaceutical and food industries. 

The translocation of a polymer chain through a 
nanochannel can be described as a three-stage process: 
the chain must first find the pore, then enter it and finally 
move through it. It is the second stage, of a polymer en- 
tering a long nanochannel, that we concentrate on here. 
Recent theoretical progress [l-Q] has shown that the en- 
try of a polymer chain into a narrow channel driven by a 
fluid flow can be regarded as a tunnelling phenomenon, 
in which the entropic cost of squeezing the chain into 
the pore is opposed by the energy gain provided by the 
driving hydrodynamic force. The outcome of such com- 
petition is a free energy barrier, which the polymer has 
to surpass in order that translocation takes place. 

According to the de Gennes blob model for conflned 
chains [5|, the barrier is overcome once the chain has 
been pushed a distance y* into the channel, at which 
point, the hydrodynamic force wins over the entropic 
pressure. Therefore, the strength and position of the 
the free energy barrier can be controlled by varying the 
driving volumetric flux, J. Increasing J has the effect of 
shifting the position of the barrier closer to the channel 
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entrance, up to a critical flux, Jc, at which the barrier 
height becomes comparable to the thermal energy and 
translocation takes place. This gives the scaling relation 



where ks is Boltzmann's constant, and T and rj are the 
temperature and the viscosity of the fluid medium. 

Remarkably, the threshold flux given by eq. ([l} is in- 
dependent of the degree of polymerisation of the chain, 
N . Physically, this occurs because the free energy bar- 
rier is overcome when a major portion of the chain is still 
outside the pore. This behaviour thus belongs to a long- 
chain regime, where the polymer always reaches the posi- 
tion of the barrier before being completely brought into 
the channel. Conversely, there is a short-chain regime, 
in which the whole chain is pushed inside the pore with- 
out reaching y* . In this case, the barrier still arises from 
the competition of the entropic pressure and the hydro- 
dynamic force as before, but with the crucial difference 
that its magnitude and position are determined by the 
length of the conflned chain. The practical consequence 
is that the translocation process becomes A^-dependent, 
a feature of potential interest, for example, in sorting 
chains according to their number of monomers. 

In this paper we shall focus on the existence of two 
regimes of translocation between long and short chains, 
and will examine the effect of the size of the chain on 
the threshold flux in each regime. Applying the energy 
barrier approach of Ref. Q we obtain a different scaling 
relation for the threshold translocation flux for long and 
short chains and test our predictions against numerical 
simulations. 

The numerical modelling of polymer-solvent dynam- 
ics has made considerable progress during the last few 
years. In particular, coarse-grained models for the fluid 
that offer a physical coupling with polymer chains have 
been developed @, 0, leading to a reliable representa- 
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tion of the dynamics of the system whilst retaining an 
efficient numerical performance. Given that the large- 
N limit relevant to the blob model is difficult to access 
directly, we shall adopt a coarse-grained representation 
introduced by Diinweg et al. [1] and described by Usta 
et al. 0, which is based on a lattice-Boltzmann model 
for the fluid coupled to a bead-spring representation of 
the polymer chain through an effective Stokes drag. This 
approach has been used to study the lateral migration of 
chains in channel flows under external forces and pres- 
sure gradients [lH- More recently, Markesteijn et 
al. used the same model to study the forced translo- 
cation of long chains into narrow pores in the limit where 
the bead-spring model is expected to give the blob-model 
behaviour. As expected, they confirmed the scaling of 
eq. (H} showing that the coarse-grained model is indeed 
able to capture the main mechanisms at play in the 
translocation process. 

The rest of this paper is organised as follows: in Sec- 
tion In] we present scaling arguments for the position 
and height of the energy barrier, and for the thresh- 
old translocation flux in the long-chain and short-chain 
regimes. In Section lllll wc describe the lattice-Boltzmann 
algorithm and the bead-spring model for the polymer 
chain, and list the set of parameters used to carry out 
the numerical simulations. Our numerical results are 
presented in Section IIVI After describing the simula- 
tion setup in Section HVAl in Section HVBI we study the 
translocation process for long chains, demonstrating the 
presence of the energy barrier. We then focus on the 
cross-over to the short-chain regime in Section HV C[ and 
on the short-chain translocation process which gives rise 
to a length-dependent threshold flux, in Section IIVDI 
Finally, in Section |V] we present the discussion and con- 
clusions of this work. 



II. ENTRY OF A POLYMER CHAIN INTO A 
NANOPORE: BLOB MODEL 

In this section we present scaling arguments to pre- 
dict the threshold flux allowing the translocation of linear 
polymer chains through a narrow channel. As we have 
anticipated above, wc will examine two regimes for the 
translocation corresponding to long and short chains. 
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FIG. 1: Schematic representation of a long blob-like polymer 
chain under confinement. The chain is confined up to a dis- 
tance y, adopting a configuration of M stacked blobs of size 
5. For a linear chain, the size of a blob scales as ^ ~ D. When 
subject to a driving flux, J, the competition between the hy- 
drodynamic drag and the entropic pressure sets a barrier to 
translocation, located at y*. Once it has reached the position 
of the barrier the polymer is able to translocate. 



chain, L, follows from the minimisation of the Flory en- 
ergy 



and reads, 
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(2) 



(3) 



In the de Gennes' blob picture [1^, the confined chain 
accommodates itself into M blobs of uniform size ^. The 
number of blobs then obeys the relation M = LD'^ /^^ . 
Within each blob the effect of confinement is unimpor- 
tant. Hence, the size of the blob scales as ^ ~ aP^^^, 
where P is the number of monomers in each blob. As- 
suming that the blobs are stacked forming a mesh, one 
can write, for the volume fraction. 
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from which it follows that for a linear chain, the size of 
the blob is comparable to the pore size, i.e., 



(4) 



A. Confined chains in equilibrium 

We start by reviewing the scaling argument of Sakaue 
et al. @ for long, blob-like linear chains. Consider a lin- 
ear polymer chain of ideal radius Rq ~ aN^/'^ , composed 
of N monomers of size a. In equilibrium, the Flory radius 
of the free chain is i? ~ aN'^/^ . 

The size of the chain changes when it is confined in 
a channel of width D < R. As depicted in Fig. [l] the 
chain stretches due to the effect of confinement. For a 
completely confined polymer, the equilibrium length of 



B. Long-chain translocation 

Wc now consider a partly confined chain in the pres- 
ence of a driving fiow. For a weak driving flux, the confor- 
mation of the chain outside the channel is close to equi- 
librium. Hence, the translocation process is controlled by 
the forces acting on the conflned part of the chain. This 
corresponds to the inside approach in the terminology of 
Sakaue et al. and is valid as long as the length of the 
pore, Lp, is larger than its thickness, D p^ . 
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The confinement of the chain has an entropic penalty 
of the order of ksT per blob [l^. Thus, for a chain 
composed of M blobs which penetrates a distance y into 
the channel, as shown in Fig. [1] there is an energy cost 



Fs ~ keT^^ - AkBT (-| 



(5) 



where we have used eq. ^ , and introduced the numerical 
prefactor A to drop the ~ symbol. As expected, the 
entropic cost increases with y, given that a larger number 
of monomers are pushed into the channel. 

Countering the entropic cost is the fluid flow, which 
tends to drag the chain further into the channel. The hy- 
drodynamic drag per blob scales as jyu^, where u = J/D^ 
is the typical velocity in the translocation direction in- 
side the channel. For M blobs, the change in free energy 
corresponds to the work done by the fluid to displace the 
chain up to a distance equal to y, 



Fh ^ -vu / M{y')^{y')dy' 
Jo 



B-qJ 



(6) 



which also increases with the penetration length, both 
because the dragging fluid has a larger number of blobs 
on which to act, and because each blob has moved further 
into the channel. As before, wc introduce a proportion- 
ality constant, B. 

Adding the contributions given by cqs. ([5]) and ([6]), it 
is possible to write the free energy change due to confine- 
ment under flow. 



{y) = Fs + Fh 



Bi]J f y 
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The competition between the entropic and hydrodynamic 
terms gives an energy barrier, AF* , located at y = y*. 
Differentiation of eq. ^ gives the position of the barrier. 



y =[t; 



(8) 



which can be substituted back into eq. ([7]) to obtain its 
magnitude. 
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where C = [A/BfB. 

In order that translocation proceeds, the chain must 
overcome this barrier. To calculate the threshold flux we 
consider the probability of translocation of the chain. 
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where kjj ~ kBT/rjD^ is a characteristic frequency and 
T,n is the observation time of the process. We define the 
threshold flux by 



P{Jc) = Pc, 



(11) 



where Pc is an arbitrary threshold probability. Inverting 
eq. (|lip leads directly to the scaling relation ([T|), 



Jc 



ksT 



As wc had anticipated, Jc does not depend on N. This 
is because the polymer reaches the position of the barrier, 
?/*, before it is completely confined in the channel, i.e., 
the translocation is triggered when A^* < N monomers 
have been pushed into the channel. To obtain the critical 
number of monomers, N*, it suffices to replace L by y* 
in eq. ([3]), from which 
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where we have used eq. We note that in Ref. @ 

the scaling for the threshold fiux was obtained by setting 
AF* ~ hgT in eq. This corresponds to a free en- 
ergy barrier whose height is overcome by pushing only 
one blob into the pore, with a corresponding number of 
monomers No = N*{y* = D) ~ {D/af^. 



C. Short-chain translocation 

The short-chain regime corresponds to the limit L < 
y* , or cquivalcntly, to A^ < A*, where the chain is 
completely confined before overcoming the energy bar- 
rier given by eq. ([9]). From eqs. (|8]) and (fT2)) . for a given 
chain length, this corresponds to weak fluxes and/or wide 
pores. Instead of being located at y*, the energy barrier 
corresponds to a penetration y — L, given that there is 
no extra cost for pushing the chain into the channel any 
further than its own length. In terms of the number of 
beads, the free energy barrier obeys AF* = AF(L(A^)), 
and follows from eq. ([7]), 



AF* 



2A* ' 



(13) 



where we have used eqs. and ([T^ to write the depen- 
dence on A^ and A^* explicitly. According to this result, 
for a given constant value of the flux, the barrier can be 
reached more easily by decreasing the number of beads. 
Conversely, increasing N has the effect of increasing the 
strength of the barrier, up to Af = A^*, where one crosses 
over to the long-chain regime, and eq. reduces to 
eq. ([9]) as expected. 

The threshold flux can be calculated from the prob- 
ability of translocation, which follows after combining 
cqs. (fTO)) and (|T3|) . Using the criterion given by eq. (|lip . 
we obtain 
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The scahng should be valid in the range No < 
N < N*, which corresponds to D < R and L < y*. 
In this range, the threshold flux increases monotoni- 
cally with the number of monomers. This behaviour can 
be traced back to the A^-dependent terms in eq. ([TS]). 
where the linear term, corresponding to the entropic cost, 
grows faster than the hydrodynamic gain, provided that 
N/N* < 1. 



III. NUMERICAL METHOD 

To further study the cross-over from long- to short- 
chain translocation discussed in Section |lTl we use the 
numerical method introduced in Ref. @ . This is a hybrid 
scheme that couples a lattice-Boltzmann fluid [H, to 
a bead-spring model for the polymer chain. Our work is 
based on the original implementation of the code, susp3d, 
which was kindly provided by the developers. Here we 
present the main features of the numerical method. For 
a more detailed description, the reader is referred to the 



original papers 
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A. Lattice Boltzmann Method 

In the lattice-Boltzmann algorithm the fluid dynamics 
follows from the evolution of the particle velocity dis- 
tribution function /j, which is defined for a discretised 
velocity set {ci}. For a given velocity vector, fi is pro- 
portional to the average number of particles moving in 
the direction of c^. 

The dynamics of /,; is given by the lattice-Boltzmann 
equation. 



Mf+ c,At, t + At)- h{r, t) = A,(r, t) + F^'- 



(15) 



which is defined over a lattice composed of nodes that 
are joined by the link vectors CiAt, where At is the 
time step of the algorithm. Here we use the D3Q19 
model, which consists of a cubic lattice with a set of 
nineteen velocity vectors in three dimensions. The lat- 
tice spacing. Ax, is uniform along the lattice axes. The 
model has three possible magnitudes of the velocity vec- 
tors, {\ci\} = (0, 1, •\/2)Ax/At. Accordingly, the fi have 
a corresponding weight that satisfies the condition 
= 1 . A suitable choice of the weights for the 
lattice model used here is ao = 1/3, ai = 1/18 and 
aV2 = 1/36. 

The dynamics expressed by eq. (jlSp is composed of 
two steps. First, the distribution function undergoes a 
collision step, where fluid particles exchange momentum 
according to the collision operator, Aj, and are driven 
by the term Ff'^*, which plays the role of a body force. 
Following this collision stage, the fi are propagated to 
neighbouring sites in a streaming step, corresponding to 
the left-hand side of the equation. 



The mapping between the lattice-Boltzmann scheme 
and the hydrodynamic equations follows from the defi- 
nition of the hydrodynamic variables as moments of the 
The local density of mass, p, momentum pt/, and the 
momentum flux, H, are given by 



fi = P, 



fiCi = pv, and ^ /^i 



H, 



(16) 

respectively. 

In the absence of external forces, the system relaxes to- 
wards equilibrium through the collision stage in eq. (|15p . 
For the lattice-Boltzmann model used in this paper, this 
is done by defining the post-collision distribution func- 
tion, /* = /; -t- Ai , which can be expressed as an expan- 
sion in the fluid velocity: 



pv ■ Ci [pvv + n"°q^*) : (c* c* -cll) 



f: = acAp + 



(17) 

where = 3 ^'"^Ax/At is the speed of sound in the 
model. Similarly, the forcing term, i^™', can be ex- 
panded in powers of the fluid velocity (cf. Ref. [31 )■ 

The post-collisional momentum flux tensor, H"°'^'*, de- 
scribes the relaxation towards the equilibrium momen- 
tum flux tensor, H""^, according to 

jjncq,* ^^^^ ^^Qncq + 1(1 + A^)(H"<='5 : 1)1, (18) 
o 



where IV""^ H - H^i, H 



cq 



+ pvv, and the bar 



indicates the traceless part of H"'''^. 

The parameters A and A^ characterise the relaxation 
timescales of the lattice-Boltzmann fluid. By performing 
a Chapman-Enskog expansion of eq. (jisp . corresponding 
to the limit of long length and timescales compared to 
the lattice spacing and the relaxation timescale of the 
fluid, the lattice-Boltzmann scheme leads to the continu- 
ity equation for the fluid density and the Navier-Stokes 
equations with second order corrections in the velocity 
for the fluid momentum. In this limit the parameters A 
and Xu map to the shear and bulk viscosities of the fluid. 



r; = -pciAt [\ + ]- 



and 



2pc2 /I 1 
4"^Ha + 2 



(19) 



(20) 



respectively. 

Solid boundaries in the simulation box are imple- 
mented using the well-known bounce-back rules |14|. 
These correspond to a reflection of any distribution func- 
tion propagating to a solid node back to the fluid node 
it came from at the streaming stage in eq. ([T5|) . As a 
consequence, a stick condition for the velocity is recov- 
ered approximately halfway between the fluid node and 
the solid node. 
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B. Polymer chain 

The linear polymer is modelled by a chain composed 
of N beads joined by freely rotating bonds. The bonds 
behave like Hookean springs, with an clastic potential 



Uci{r) = k{r 



(21) 



where r is the separation between adjacent beads, a is 
the equilibrium length of the bond and k is the elastic 
constant of the spring. 

The short-range excludcd-volume interactions are 
modelled by a truncated DLVO potential, 



C^DLVO = Uq 



exp(-KDH?') 



(22) 



where kdh is the Debye-Hiickel screening length and Uq 
is an amplitude. 

The position vector of the i-th polymer bead, Xi, 
evolves in time according to 



mxi 



:U + Fi. 



(23) 



where the summation term includes the excluded- volume 
interactions with all beads and the elastic interactions 
with neighbouring beads. The term Fi contains the vis- 
cous hydrodynamic force that couples the bead to the 
lattice-Boltzmann fluid. 



Fi = -^oixi - v{xi)) + Fl. 



(24) 



This expression includes the Stokes drag, = STrr^rn, 
with ru being the hydrodynamic radius of the beads, 
and a random force Fl that satisfies the fluctuation- 
dissipation relation. 



{Fl{t)Fl{t'))^2kBm5{t~t')\. 



(25) 



While the polymer beads move in continuous space, 
the lattice-Boltzmann fluid is only defined at the lattice 
nodes. In order to calculate the coupling force, Fi, the 
model uses a linear interpolation scheme to estimate the 
fluid velocity, i/, at the position of the bead, Xi. The 
discretisation of the lattice gives an effective hydrody- 
namic radius, r^, which differs from the input hydro- 
dynamic radius, rjj. Here we follow the same procedure 
as in Ref. Q, and choose rn in order to obtain the de- 
sired value of . Once the hydrodynamic force has been 
exerted on the monomer, momentum conservation is en- 
forced by exerting a force of equal magnitude back onto 
the fluid. 



C. Parameter Values 

Our objetive is to carry out numerical simulations of 
polymer chains in the limit where the bead-spring model 
gives the blob-model behaviour presented in Section |lTl 




FIG. 2: Schematic representation of the simulation box. 



Such regime has been validated when the monomer size 
matches the blob size, a ~ ^, for linear chains, where 
^ ~ D (3]. Therefore, the bead-spring model should give 
the blob-limit behaviour for a ~ D. On the other hand, 
the hydrodynamic radius of the beads, , is limited to 
small values, corresponding to the point-particle coupling 
between the chain and the lattice-Boltzmann fluid Q. 
With these conditions in mind, we fix the bead diameter 
to 2r^ = Aa;/2, while a = Ax and D = 2Ax-. 

The remaining model parameters are chosen as follows: 
we work at a fixed temperature, ksT = 0.1. To pre- 
vent chain crossings the spring constant is taken as fc = 
SOOfesT/Ax^. Parameter values for the excluded- volume 
potential arc fixed to — ksTAx and kdh = 80/Aa;, 
which ensure that inter-monomer repulsions are larger 
than k bT a.t distances comparable to . To match the 
hydrodynamic diameter to the effective bead diameter, 
wc fix ru = 0.32Aa; following the calibration procedure 
presented in Ref. |9|. 

In order to resolve the hydrodynamics correctly, one 
needs to ensure that the distribution function relaxes on 
a faster timcscalc than the diffusive timcscalc of the poly- 
mer chain. This condition can be satisfied by setting the 
parameters in the collision operator in eq. (|15p to A = — 1 
and Xi, = —1. We set the fluid density to p = 36, and 
the timestep and lattice spacing in the lattice-Boltzmann 
fluid to At = 1 and Aa; = 1. Using these values, the dy- 
namic viscosity in simulation units is 77 = 6. 



IV. NUMERICAL RESULTS 
A. System Geometry and Initial Conditions 

The geometry of the system is depicted in Fig. [5J We 
consider two identical rectangular ducts of dimensions 
= 108, Ly ~ 14 and = 10, separated by a square 
pore of width D = 2 and length Lp = 24. Solid walls 
are imposed in the y and z directions, while periodic 
boundary conditions are enforced in the x direction. 

For small Reynolds numbers, one expects that the vol- 
umetric flux, J, caused by applying a uniform body force 
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FIG. 3: Volumetric flux, J, as a function of the body force, 
showing that Darcy's law holds. 

to the fluid, /, obeys Darcy's Law: 

J=-L (26) 
V 

where S and k are the local cross-sectional area and per- 
meability of the medium and rj is the dynamic viscosity 
of the fluid. 

In order to verify the validity of eq. (j26p . we measure 
the local flux for different applied forcings. The results, 
illustrated in Fig. [31 show a linear dependence of J on 
/ as expected. A linear fit to the data gives kS ~ 238. 
For a given forcing, we find that the flux does not vary 
significantly whether it is measured inside or outside the 
pore thus confirming that, for the range of forcings con- 
sidered here, the lattice-Boltzmann fluid behaves as an 
incompressible liquid. Given that we are interested in 
the scaling of the critical flux, we use values of the forc- 
ing that give a ratio between hydrodynamic and thermal 
effects in the range < ryJ/fc^T < 1.5. 

We initially tether the polymer chain at a position 
y = yo inside the channel as depicted in the inset of Fig.SJ 
and let the system equilibrate for 5x 10^ time steps. The 
equilibration period allows the relaxation of the chain. 
Evidence for this is presented in the same figure, which 
shows the number of beads inside the pore, Ny^, as a 
function of yo for = 128. For a confined chain in equi- 
librium, we have, from eq. ([3|), Ny^ ~ (yo/D){D/a)^/^ . 
Our results show that the number of beads increases lin- 
early with the tethering position, as expected. 

B. Energy Barrier 

In order to confirm that the polymer translocation 
is controlled by surpassing an energy barrier, we first 
perform simulations of long polymer chains subject to 
a fixed driving flux while varying the initial tethering 
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FIG. 4: Number of beads in the confined part of the chain as 
a function of the initial tethering position. Inset; schematic 
representation of the tethered chain. 



position, yo- We consider different values of yo in the 
range < yo < Lp, and set the number of monomers to 
— 128, thereby ensuring that the major part of the 
chain lies outside the pore. 

Simulations are carried out by letting the chain equi- 
librate as before. Once the chain has equilibrated, it is 
released from its tethering point and the body force is 
applied uniformly to the fluid. The chain then is either 
carried down the pore by the underlying fluid flow and 
eventually translocates to the opposite duct, or is ejected 
from the pore back into the original chamber. Simula- 
tions are run for 5 x 10^ timesteps, which is a sufficiently 
long timescale to identify successful or failed transloca- 
tion events. 

Fig. \5\ shows the probability of translocation of the 
chain as a function of yo for three different values of the 
imposed flux. We observe a clear transition from non- 
translocating (P = 0) to translocating (P = 1) chains as 
yo is increased. This confirms the presence of an energy 
barrier, which is progressively approached as the chain is 
pushed further into the channel. The probability curves 
are shifted to the left as one increases the imposed flux, 
J. This indicates a shift of the position of the barrier, 
y* , closer to the pore entrance caused by a higher driving 
hydrodynamic force. 

As a criterion we define the position of the barrier as 
P{y*) ~ 1, and plot the measured values of y* as a func- 
tion of ksT/riJ in the inset of Fig. [5] As expected from 
eq. ([8]), we observe a linear growth. A fit of the data 
shown in the figure to the function y*{x) ~ {A/B)Dx 
gives an estimate of the numerical prefactor in eq. ([5]), 
A/B ~ 4.4. 

The imposed flux sets the position and height of the en- 
ergy barrier. Given that we fix the chain to a prescribed 
position inside the channel, the height of the barrier is 
reduced by an amount determined by yo- This can be 
easily included in the analytical model by replacing y by 
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FIG. 5: Probability of translocation as a function of yo- In- 
set: linear growth of the position of the energy barrier with 
keT/riJ. 



y — yo in eq. ^ and performing the integration in eq. ^ 
from y' = yo io y' = y. The resulting free energy barrier 
reads 



^F* A fy*\ ( fyo 



(27) 



where g{x) = —2x+x^. From this expression and eq. PU)) 
we can calculate the probability of translocation. This 
predicts logP ^ ls.F* /ksT as observed in Fig. El As 
expected, all data points collapse onto the same curve, 
which is linear in the free energy. A fit to the data 
gives an estimate for the numerical prcfactor in eq. (|27p . 
A~ 1.8, which together with the measured value of A/B, 
gives the amplitude of the hydrodynamic contribution to 
the free energy in eq. dH]), i? — 0.41. 

The inset in Fig. [5] shows how the theoretical predic- 
tion accurately captures the main features of the translo- 
cation process. The probability of translocation increases 
to unity as yo/y* — >■ 1. The spread of the curves is dic- 
tated by the ^ ksT/riJ prefactor in eq. (P7)) . increasing 
the likelihood of translocation at a given yo/y* for larger 
fluxes, or lower temperatures. 
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FIG. 6: Logarithm of the translocation probability as a func- 
tion of the free energy barrier magnitude. Symbols correspond 
to numerical simulations. The solid line is a linear fit to the 
data. Inset: translocation probability as a function of yo/y* . 
Numerical results (symbols) show a good agreement with the 
theoretical prediction (solid lines). 



1' 
0.8 
0.6 
0.4 
0.2 











♦ 


— <^ 




















— ^ 








1 ^''^ i 


— ■ — 


■ 



16 



24 



32 



40 
N 



48 



56 



64 



nJ/l<BT=0.26 h 



nJ/l<BT=0.64 h 



nJ/l<BT=0.43 I— •— I nJ/l<BT=0.86 h- ♦— I 
nJ/l<BT=0.50 h-A^ 



FIG. 7: Probability of translocation as a function of the num- 
ber of beads in the chain at different values of the imposed 
flux. Simulation results correspond to symbols, while the the- 
oretical model is indicated by the solid lines where we set 
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C. Cross-over from long to short chains 

So far we have discussed the translocation of long 
chains, which can reach the position of the barrier while a 
significant amount of monomers remain outside the pore. 
The threshold flux is thus controlled only by the number 
of monomers that it takes to reach the barrier, N* , and 
not by the total number of monomers in the chain, N . 

This picture breaks down for shorter chains, where our 
theory predicts that the scaling of the critical flux be- 
comes dependent on N . In order to explore this effect, 
we have carried out simulations at a fixed initial tethering 



position, j/Oj a-nd driving flux, while decreasing the total 
number of beads in the chain, N . We fix the initial teth- 
ering position to yo = Lp/2 = 12, for which the number 
of beads inside the pore after equilibration is A^y^ ~ 18. 
We carry out the simulations in the range 16 < A^ < 64, 
where we expect to observe the cross-over. 

Figure [7] shows the probability of translocation as a 
function of A^ for five different values of the driving flux. 
The probability always increases with the applied flux as 
expected. For large A^ the probability is insensitive to 
the chain length, over a range that persists to smaller A^ 
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FIG. 8; Probability of translocation as a function of the im- 
posed flux for different number of monomers in the chain. 



as the flux increases. For smafler N the probabihty of 
transloeation increases strongly with decreasing N until 
it saturates at ~ 20 to P — >■ 1, where all polymer 
chains translocate even under the weakest flow. 

This behaviour is in agreement with our theoretical 
model. For large N, the critical number of beads to 
overcome the barrier follows from eq. (|12p . and obeys 
N* ~ {D/a)^/'^{kBT/T]J). Therefore, the cross-over to 
the short-chain regime [N < N*) is observed at larger 
chain lengths as J is decreased. From eq. (fT2|) we also 
have N* ~ y*, from which it is possible to estimate the 
critical number of beads corresponding to each driving 
flux by interpolating the data shown in Fig. [H Once 
each value of N* is known, we can determine the energy 
barrier as a function of N from eq. (|13p . In order to in- 
clude the effect of yo, we follow the same procedure as 
that used to obtain eq. (P7)) . This gives 



AF* _ ^( ^Y^^ 
ksT ~ \d) 



N 



N 



1-i 
2N 



-N* 
2 



N* 



(28) 

Using this expression for the free energy barrier an es- 
timate for the probability of translocation follows from 
eq. (Uni). The observed saturation length ~ 20 is in- 
deed in good agreement with our estimate Ny^ ~ 18, at 
which the free energy barrier given by eq. (j28|) vanishes. 
As shown in Fig. [71 we obtain a good agreement between 
the theory and simulations. 



which the probability curves cross over from complete 
rejection of the chains to complete translocation. 

Figure [8] shows the resulting probability curves. For 
> 32 the iV-independence of the long-chain regime is 
recovered, as curves fall on top of each other. Our results 
in the long-chain regime are consistent with recent exper- 
imental studies of polymer translocation in nanochannels 
carried out by Beguin et at. [l^. In their experiments 
they measure the rejection coefficient of chain transloca- 
tion, which is related to the probability of translocation 
by 7?. ~ 1 — P. Their experiments give the same smooth 
transition from rejection to translocation of the chains in 
the range < rjj /ksT < 1, very close to our simulation 
results. For iV < 32 the probability curves shift upwards 
and become increasingly plateau-like as the number of 
monomers is decreased, indicating the cross-over to the 
short-chain regime. 

In order to quantify the cross-over from the short- 
to the long-chain regime, we use the criterion given by 
cq. (fTT|) to interpolate Jc from each of the curves shown in 
Fig.[S] The threshold flux is shown in Fig.lHlas a function 
of N for six different values of the threshold probability. 
As expected, the overall behaviour does not depend on 
the particular choice of Pc . The threshold flux saturates 
for large N, corresponding to the long-chain regime, and 
shows a marked decrease as N Ny^, where the chain 
becomes completely confined inside the channel. 

Using eqs. (jlOp . (fTTj) and (|28l). we can estimate the 
threshold fiux as a function of the number of monomers 
in the chain. 



■ilJc 



_2^ 



10/3 



(29) 



[^f'iN-Ny„) + 



log(Kprni) ' 

A 



Both prefactors, A and B, depend on the intrinsic 
properties of the chain and not on the particular translo- 
cation regime. We therefore set them to the values that 
we have measured previously. 

Figure [9] shows the resulting comparison between the- 
ory and simulation. Remarkably, we find a good quanti- 
tative agreement by using no new fit parameters in our 
theoretical prediction. 



V. SUMMARY AND DISCUSSION 



D. Critical flux for small chains 

We now turn our attention to the dependence of the 
critical flux, Jc, on the number of monomers in the short- 
chain regime. To find the threshold flux, we carry out 
simulations of the translocation process at a fixed number 
of beads while increasing the flux. We fix ?/o = 12 as 
before, and consider the range < rjJ/kBT < 1, for 



The results presented in this paper show that the 
translocation of polymer chains through narrow pores 
exhibits two different regimes depending on the length 
of the chain. As previously proposed in Ref. [H, our 
numerical results show that the translocation process is 
controlled by overcoming a free energy barrier charac- 
terised by a critical penetration of the polymer inside 
the pore, y*. For long chains, this length scale always 
outruns the polymer length under confinement, L. Thus, 
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FIG. 9: Threshold flux as a function of the number of 
monomers in the chain. Pc is the choice of threshold prob- 
ability. Simulation results: symbols. Theoretical prediction: 
solid lines. 



the translocation process is controlled by y* , regardless of 
the number of monomers in the chain, A^. Conversely, for 
small chains, where the confined length of the polymer 
is smaller than y* , we have shown that the transloca- 
tion process is controlled by L, and the process becomes 
A'^-dependent. A qualitatively similar conclusion should 
apply to the passage of the chains through two dimen- 
sional narrow slits |16| . 

For even smaller chains, where L is comparable to the 
size of the pore, D, and for very small fluxes {qJ/kBT ~ 
10~^), we have observed a reduction of the probability 
of translocation of the chains. This is caused by the ef- 
fect of diffusion of the chain inside the pore, which can 
favour the ejection of small chains at very low forcings. 
However, this regime falls out of the free energy barrier 
picture presented in this paper, and we therefore leave it 
for future exploration. 

An application of the dependence of the translocation 
probability of short chains on N is the potential sorting 
of the chains according to their number of monomers. In 
the short-chain regime the polymer must be completely 
pushed in before translocating through the pore. Given 
that only a small number of blobs can be pushed into the 
pore by fluctuations, one expects that only very short 
chains translocate by virtue of thermal effects in this 
regime. For longer chains, but still in the short-chain 
regime, the translocation must be assisted, for exam- 
ple, by a molecular motor that is able to push the chain 
deeper into the pore. 

In an aqueous solution at room temperature, the crit- 
ical flux in the long chain regime follows the scaling of 



eq. Il]) and is of order Jc ^ 10~^® m^ s~^. For a microme- 
tre sized channel the corresponding velocity, Vc — Jd , 
is of order 10~^ m s^^. Such a small value suggests that 
chains can easily translocate through microfluidic cham- 
bers at normal operation velocities, which are typically 
of centimetres to metres per second. However, recent 
experimental measurements of cytoplasmic streaming in 
micrometre cell channels |17| show that the streaming ve- 
locities are in the range of 10 /xm s~^. For weaker fluxes 
than Jc, one enters the short-chain translocation regime. 
Therefore, it is feasible that the translocation regime pro- 
posed in this paper might be at play in biological systems 
as a regulator of selective translocation. 

Here we have considered polymers in a good solvent 
using the DLVO potential to implement short-ranged 
excluded- volume interactions with the wall. Understand- 
ing the details of electrostatic effects, with regards both 
to the interaction with the wall and the interactions be- 
tween segments, and how they are affected by the salt 
concentration, is an important issue for future work, rel- 
evant, in particular, to the dynamics of DNA in practical 
applications. 

Finally, we comment on the feasibility of an experi- 



mental confirmation of our prediction. Bcguin et at. [15 1 
have recently performed an experimental study of the 
flow-driven translocation of hydrosoluble polymers into 
nanopores. In their experiments, they consider chains 
whose radius of gyration is i? ~ 84 nm, and which are 
forced across pores 23 — 35 nm in radius and 6.5 ^m 
in length. Their results for the rejection coefficient, 
72. = 1 — exp (AF*/fcBT), show a good agreement with 
the long-chain limit of the de Gennes model, where they 
use a similar expression to eq. ([9]) . According to our pre- 
diction, the short chain regime would be observable for 
chains whose length under confinement is smaller than 
the position of the barrier. This traduces into a cross- 
over radius of gyration R* ~ D{kBT /riJ)^/^ . Taking 
Z? = 70 nm and rjj /ksT w 0.8, which correspond to ex- 
perimental conditions reported in ref. (isj . we estimate 
R* K, 80 nm, which is close to the radius of gyration 
of the polymers used in the experiments. This supports 
the feasibility of future experimental work to verify the 
theoretical predictions presented in this paper. 
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